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Abstract. As part of condensed-matter physics, the field of Anderson localization concerns 
the study of conductance of electrons in a random medium. 

We consider the discrete random Schrodinger operator on the integer lattice Z 2 . Based 
on a recent mathematical result we introduce a new numerical approach to the Anderson 
(de)localization problem. As an application we show numerically the following unexpected 
result: For small disorder, this random operator allows extended states with positive proba- 
bility. This approach eliminates potential problems due to boundary effects, and the numer- 
ical part is rather simple compared with other experiments in the field. Further, we provide 
a mathematical derivation of a quantity closely related to Thouless' dimensionless scaling 
parameter. No new information is gained regarding the energy regimes at which diffusion 
occurs. 



1. Introduction 

In 1958 P.W. Anderson [2] suggested that sufficiently large impurities in a semi-conductor 
could lead to spatial localization of electrons, called Anderson localization. The field has 
grown into a rich theory and is studied by both, the physics and the mathematics community. 

We consider the discrete random Schrodinger operator on the integer lattice of dimension 
two H w = — A + Yliez 2 u i < ' > &i > $i an d assume that the random variables Wj are inde- 
pendent identically distributed with uniform distribution in [— c, c], i.e. each number in the 
interval is attained with equal probability. This operator describes the situation where the 
atoms of the crystal are located 'near' the integer lattice points Z 2 . 

Many approaches to Anderson localization are based on scaling theory. A dimensionless 
scaling parameter [T3] is the sole indicator whether or not the system exhibits Anderson 
localization. The parameter g is the ratio between the Heisenberg time (the maximum time 
that a wave packet can travel inside a finite region before revisiting the same location) and 
the Thouless time (the time it takes an extended state to arrive at the boundary of a finite 
region). Roughly speaking, for g < 1 we have localization, and otherwise the existence of 
extended states. Our results are in contradiction with much numerical and experimental 
work using the scaling theory (see e.g. [U HUJ [H]) which indicates that one should expect 
some "marginal localization" in the form of so-called quasi-extended states in dimension two 
for small disorder. 

However, mathematically rigorous proofs for the dynamical localization of the discrete 
random Schrodinger operator on / 2 (Z d ) are only known in dimension one at all strengths 
El [7] and in dimension two and higher for disorders c above a certain threshold (3J EJ [13] . 
For the discrete random Schrodinger operator on ? 2 (Z d ) there is no proof, whether or not 
diffusion occurs in dimension two and higher for small disorder c. Many results are known 
for related operators, e.g. the discrete random Schrodinger operator on the tree and other 
structures. 

The author is partially supported by the NSF grant DMS-1101477. 

1 



2 



CONSTANZE LIAW 



Based on a recent mathematical result PQ - which is based on the study of rank one 
perturbations - the author designed a new numerical approach to indicate conductance in 
random media. The key to making this method numerically feasible is that it suffices to track 
the evolution of just one vector under the repeated application of the random operator. The 
details of this numerical experiment, the mathematical background, and related information 
can be found in |12j . 

In Section [2j we explain the idea behind the experiment and how it indicates that the 
discrete random Schrodinger operator in two dimensions does in fact exhibit diffusion for 
small disorders c < 0.7 with positive probability. Via variations of the code we study the 
evolution over time of the energy of a wave packet initially located at the origin, as well 
as the dependency with respect to c of the location of the energy of such a wave packet in 
Section [3l We show how those mathematical and numerical results are related to the Thouless 
parameter, and conclude the existence of extended states in dimension two for small disorder. 

The method within can be applied to many operators from a larger class of discrete Ander- 
son models, so-called Anderson-type Hamiltonians [9j. The author expects a similar result for 
the discrete random Schrodinger operator in dimension three and higher but has not investi- 
gated those cases. However, this method cannot not (yet) be applied to two other important 
random Hamiltonians: Due to memory restrictions, the method is not suited well for the 
discrete random Schrodinger operator on the tree. And the mathematical result cannot be 
applied for the Rademacher potential where the random variables attain the values ±1 with 
equal probability. 



2. Key steps of the new approach 

Consider the discrete random Schrodinger operator H w on Z 2 (Z 2 ) with identically dis- 
tributed random variables Uij with uniform distribution in [— c, c]. Fix the vectors <5oo £ Z 2 (Z 2 ) 
and 5u 6 Z 2 (Z 2 ), namely 



(2.1) 



Joo 



. . 











.. 


. . 











.. 


. . 





1 





.. 


. . 











.. 


. . 











.. 



'11 



.. 











. . 


.. 











. . 


.. 











. . 


.. 








1 


. . 


.. 











. . 



Consider the distance D™ c between the unit vector 6n and the subspace obtained by taking 
the span of the vectors {Soo, H w 5qq, H^Sqo, • • • , H^Sqq}- In other words, let 

D" c := dist^i, span{il^oo : k = 0, 1, 2, . . . , n}). 

Deep mathematical results (see Corollary 3.3 of [12]) imply that dynamical derealization 
takes place, if we can find a disorder c > for which the distance D™ c does not tend to zero 
as n approaches oo, i.e. -D WjC > where 



D, 



— lim D n 

n— >oo ' 



It is worth mentioning that this mathematical result is a 'one-sided' implication. In par- 
ticular, if D u c = 0, then we cannot conclude localization. 
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Figure 1. As a function of c we show y WtC (x's; larger function values) and 
L WjC (circles; smaller function values). We conclude D WjC rj y Wfi > L WjC > 
for c < 0.7. 

Some fairly simple arithmetic shows that D wc is given by the following expression 
(2.2) D u , c 
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where {mo, mi, . . . , m n } form the orthogonal basis of the linear subspace spanned by H^Sqq 
for k = 0, 1, 2, . . . , n of Z 2 (Z 2 ) obtained by the Gram-Schmidt process (without normaliza- 
tion), and || • || 2 denotes the Euclidean norm. 

In the numerical experiments, in order to exclude the possibility that the distance D™ c 
tends to zero logarithmically, we re-scaled time and applying the negative exponent to the 
horizontal axis which satisfied the least squares property for line approximation. The inter- 
section y U)C of the approximating line with the vertical axis is an estimate for the distance 
Ad,c Further we considered the lower estimate L W)C for y u>c obtained by taking the lowest 
y— intercept lines through any two consecutive points. We repeated the experiment for many 
realizations of the random variable u and chose the smallest values of y w c and L u c obtained 
from all randomizations. Further, we considered a wide range of disorders c. 

We include Figure [T] (see p2]) which shows these smallest values of j/ WjC and L W)C as a 
function of disorder c. 

3. Connection with the scaling parameter and existence of extended states 

We explain how the derealization result described in Section [2] implies the existence of 
extended states in the sense of the Thouless criterion. 
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Figure 2. The figure shows the energy distribution of m3ooo as a function of 
distance from the origin for several values of disorder c. 



The numerator < m^, 5u > in formula (2.2) is precisely the (1, 1)— entry of the vector m^, 
while the denominator is the total energy of H^5qq after successive orthogonalization (we 
only consider the new information of each state). 

Roughly speaking, the numerator is small often precisely if the Heisenberg time is large, 
i.e. it takes a long time to return to the same finite region. And if much of the energy of the 
state rrik is located away from the origin, then the Thouless time is small. 

In the case of small c, Figure [2] (see [12]) shows that most of the energy of is located 
'away' from the origin almost surely. In particular, < mfc,<5ii > is small often and is 
relatively large for all k. Indeed, a large part of the energy is almost surely located on the 
outermost diamond which the orthogonalization process leaves unchanged. In particular, the 
norm ||mfc|| remains comparably large. This implies that the Heisenberg time is large and 
the Thouless time is small, yielding a scaling parameter > 1 and the existence of extended 
states (in the physicists' sense), i.e. conductance of electrons. As was mentioned above, our 
method of computing a lower estimate of Mj,c also indicates dynamical derealization and 
hence the existence of extended states (in the mathematical sense meaning states that do not 
decay exponentially, and even absolutely continuous spectrum). 

As the disorder increases, Figure [2] shows that a larger part of the energy remains near the 
origin, so that < m^, 8n > is away from zero more often and the ratio < "! fc ' j. 1 ^ is generally 

II 2 

larger than for small c. However, in this case a large part of the energy remains near the 



origin. The Euclidean norm Hrn^^ in the denominator of (2.2) has little to do with the 
Thouless time. Therefore, the relationship between the ratio in formula (2.2) and the scaling 
parameter is not directly obvious. This is in agreement with the fact that our new method 
cannot be used to show localization. 
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